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One of the robust features found in simulations of core-collapse supernovae (SNe) is the prompt 
neutronization burst, i.e. the first ~ 25 milliseconds after bounce when the SN emits with very 
high luminosity mainly Ve neutrinos. We examine the dependence of this burst on variations in the 
input of current SN models and find that recent improvements of the electron capture rates as well 
as uncertainties in the nuclear equation of state or a variation of the progenitor mass have only 
little effect on the signature of the neutronization peak in a megaton water Cherenkov detector for 
different neutrino mixing schemes. We show that exploiting the time-structure of the neutronization 
peak allows one to identify the case of a normal mass hierarchy and large 13-mixing angle i9i3, where 
the peak is absent. The robustness of the predicted total event number in the neutronization burst 
makes a measurement of the distance to the SN feasible with a precision of about 5%, even in the 
likely case that the SN is optically obscured. 



PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

Despite the enormous progress of neutrino physics in 
the last decade, many open questions remain to be solved. 
Among them are two, the mass hierarchy — normal ver- 
sus inverted mass spectrum — and the value of the 13- 
mixing angle i?i3, where the observation of neutrinos from 
a core-collapse supernova (SN) could provide important 
clues [E m ■ The neutrino emission by a SN can be 
divided schematically into four stages: Infall phase, neu- 
tronization burst, accretion phase, and Kelvin-Helmholtz 
cooling phase. The bulk of SN neutrinos are emitted in all 
flavors during the last two phases with small differences 
between the and z^^^t- spectra j^i^i^]. Moreover, the 
absolute values of the average neutrino energies as well as 
the relative size of the luminosities during the accretion 
and cooling phases are not known with sufficient preci- 
sion. As a consequence, a straightforward extraction of 
oscillation parameters from the SN neutrino signal during 
the accretion and cooling phase seems hopeless. 

An alternative is the use of observables that do not rely 
on SN parameters. Such observables require only that the 
initial neutrino fluxes -Pf are different functions of energy 
and time, F°^^^^{E,t) ^ FS^^^jE,t), where x = {fi, r}. 
Since the interaction of neutrinos with matter depends on 
their flavor, identical energies and luminosities for Ve, i^e 
and would require a conspiracy of interaction rates and 
chemical composition inside the neutrinospheres. Exam- 
ples for such observables are the modulations in the SN 
neutrino signal caused by the passage of the neutrinos 
through the Earth or by the propagation of shock 
waves through the SN envelope [1,13 • If the mixing an- 
gle i?i3 is known to be large, sin^ -dis > 10~^, an ob- 
servation of Earth matter or shock wave effects in the 
experimentally most important De channel would imply 
a normal or inverted mass hierarchy, respectively. If how- 



ever the value of z?i3 is not known, a degeneracy exists 
between the case of a normal mass hierarchy and large 
z9i3 (scenario A, cf. table 1) and the case of small z9i3, 
sin^ i?i3 < 10~^, and any hierarchy (scenario C): scenario 
A and C both predict the same signature in a water 
Cherenkov detector. 

A different way to extract reliable information about 
neutrino mixing parameters is to use characteristics in 
the neutrino emission of SNe that are model indepen- 
dent. One of the most robust features of numerical SN 
simulations is the so-called neutronization ly^ burst |lClj| . 
which takes place during the first ~ 25 ms after the core 
bounce. The small number of events expected during this 
time period is compensated by the moderate dependence 
of the Vf, burst on physical parameters like the progen- 
itor mass or details of the SN models. In Sec. II, we 
discuss the astrophysical aspects of this burst in detail, 
emphasizing the robustness of the neutrino luminosities 
against variations in the input of the SN models. In 
Sec. Ill, we study the signature of the neutronization 
peak in a megaton water Cherenkov detector for differ- 
ent neutrino mixing schemes. We argue that exploiting 
the time-structure of the neutronization peak allows one 
to identify the neutronization burst even if the SN is not 
visible in the optical. A non-observation of the neutron- 
ization burst identifies the case of a normal mass hierar- 
chy and large 13-mixing angle 'i?i3 (case A), thus breaking 
the degeneracy between the neutrino mixing scenarios A 
and C. Moreover we find that for a given neutrino mixing 
scenario the systematic uncertainty due to unknown SN 
parameters affects only little the total number of events 
in the neutronization burst. As we discuss in Sec. IV, this 
robustness of the theoretical prediction makes a measure- 
ment of the distance to a SN located at 10 kpc feasible 
with a precision of about 5%. Such an accuracy is compa- 
rable to optical methods using the SN light curve, which 
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have an error between 5 and 10%. If the SN is optically 
obscured, measuring its distance through the v^. burst is 
crucial for estimating the total binding energy released or 
to limit the strength of the gravitational wave signal (Tl| 
emitted by the SN. Finally, we summarize our results in 
Sec. V. 



II. ASTROPHYSICAL ASPECTS OF THE 
NEUTRONIZATION BURST 

Modern supernova models with sufficiently detailed 
treatment of the neutrino physics have in common 
the existence of a "prompt burst" of electron neutri- 
nos El 111 El- This breakout pulse is launched at 
the moment when the newly formed supernova shock 
that races down the density gradient in the collapsing 
stellar core reaches densities low enough for the ini- 
tially trapped neutrinos to begin streaming faster than 
the shock propagates E3- In the shock-heated matter, 
which is still rich of electrons and completely disinte- 
grated into free neutrons and protons, a large number of 

are rapidly produced by electron captures on protons. 
They follow the shock on its way out until they are re- 
leased in a very luminous flash, the breakout burst, at 
about the moment when the shock penetrates the "neu- 
trinosphere" and the neutrinos can escape essentially un- 
hindered. As a consequence, the lepton number in the 
layer around the neutrinosphere decreases strongly and 
the matter neutronizes [13 ■ Because of the high tem- 
peratures behind the shock, electron-positron annihila- 
tion, nucleon-nucleon bremsstrahlung |l3j |. and, when 
Ve become more abundant, also neutrino-pair conversion 
VeVe — > i^/^.ri^/i.T ;5| are efficient in creating muon and 
tau neutrino-antincutrino pairs ^h^. The luminosities of 
the latter therefore begin rising steeply immediately af- 
ter shock formation. In contrast, the luminosity of Vf, 
increases more slowly. On the one hand this is due to the 
fact that the abundance of positrons and therefore the Vg, 
production by e+ captures is rather low as long as elec- 
trons are still highly degenerate, on the other hand pair 
creation of and is also suppressed by the high abun- 
dance of z^e during the burst phase and the corresponding 
fermion blocking in the Ve phase space. 

These facts can be verified from Figs. and [3 The 
rms energies shown in the latter figure are defined by 



/o de/_^ d/x/^(e,/x)e5 
\ /o°°de/^j^dAi^(e,M)e^ 



(1) 



with /i/(e, /i) being the neutrino phase space distribution, 
which is a function of the neutrino energy e and the co- 
sine, /X, of the angle of neutrino propagation relative to 
the radial direction. The results presented in the plots 
were obtained by core collapse simulations in spherical 
symmetry with the neutrino-hydrodynamics code devel- 
oped by Rampp and Janka |17| , employing a solver for the 



energy-dependent moments equations of neutrino num- 
ber, energy, and momentum and an approximative treat- 
ment of general relativity that yields good agreement 
with fully relativistic simulations, in particular during 
the collapse and early postbounce phases 1^. 

A local minimum in the Ue luminosities occurs shortly 
after the formation of the shock at core bounce {t = 0) 
and before the neutronization burst. It is caused by the 
shock first compressing matter from a semi-transparent 
state to neutrino-opaque conditions before the post- 
shock layer reexpands to become neutrino transparent 
and to release the neutronization neutrinos E3|- Per- 
forming simulations for a variety of progenitor stars be- 
tween 11.2 Af0 and 25 M0 from different stellar evolu- 
tion modelers |l9j |. we have confirmed the uniformity of 
the radiated neutrino luminosities and rms energies in 
the first 20 ms after bounce (Figs. ^ and [3 left pan- 
els) that was also seen in other recent simulations with 
neutrino transport being described by a solution of the 
Boltzmann equation or its moments equations 0, |23| ■ 
The prompt neutronization burst has a typical full width 
half maximum of 5-7 ms and a peak luminosity of 3.3- 
3.5x 10^"^ ergs^^. The striking similarity of the neutrino 
emission characteristics despite of some variability in the 
properties of the pre-coUapse cores is caused by a regu- 
lation mechanism between electron number fraction and 
target abundances (protons and nuclei) for electron cap- 
tures 21, 22], which establishes similar electron fractions 
in the inner core during collapse. This leads to a conver- 
gence of the structure of the central part (of roughly a 
solar mass) of the collapsing cores and only small differ- 
ences in the evolution of different progenitors until shock 
breakout Differences of the core size and of the 

density profile in the outer part of the iron core and be- 
yond lead to different mass infall rates at late times when 
the shock has reached neutrino-transparent layers. This 
implies different accretion luminosities and thus causes 
a progenitor-dependent strong variation of the neutrino 
emission characteristics after the v^. luminosity has lev- 
elled off from the prompt burst. 

Only recently improvements in the treatment of elec- 
tron capture rates on nuclei during the late phases of 
stellar evolution and core collapse have become available 
which remove shortcomings of the widely used indepen- 
dent particle model in which electron captures are sup- 
pressed by Pauli blocking for nuclei with TV > 40 |2ll| . 
This typically happens at a density of some lO^'^gcm"'^ 
above which electron captures on free protons govern the 
evolution of the electron fraction. The improved rates 
for core collapse are based on Shell Model Monte Carlo 
calculations of nuclear properties at finite temperatures, 
complemented with a Random Phase Approximation for 
the electron capture rates of a wide sample of nuclei in 
the mass range between A = 65 and A—l\2 with abun- 
dances given by nuclear statistical equilibrium |23j . In 
supernova simulations with these improved rates electron 
captures by nuclei dominate over capture on free protons, 
and interesting changes were found during core collapse. 
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FIG. 1: Luminosities as functions of time for Ve (top), Ve. (middle) and heavy-lepton neutrinos (bottom). In the left column 
results for different progenitor stars between 11.2 M© and 25 (left column; the progenitor mass is indicated by the number 
after the first letter of the model name ITg|l 'l are shown, in the middle column for simulations with the new treatment of electron 
captures by nuclei during stellar core collapse according to Langanke, Martmez-Pinedo and Sampaio (LMS; red solid and dotted 
lines) compared to the traditional description (black lines) in case of a 15 Mq and a 25 Mq star. The right column shows results 
for three different nuclear equations of state applied to the collapse of a VoMq progenitor (see text for more details). The 
luminosities are given for an observer at rest, evaluated at a radius of 400 km with a corresponding time retardation of about 
1 ms. Time is normalized to the moment of shock formation defined by the instant when the entropy behind the shock first 
exceeds a value of 3 fcs per nucleon. 



bounce, and postbounce evolution j23ll24j. In the panels 
of the middle columns of Figs. ^ and |3 one can see the 
corresponding differences in the neutrino emission prop- 
erties for simulations of a 15 Mq and a 25 Mq progeni- 
tor with the new capture rates accor ding to Langanke, 
Martmez-Pinedo and Sampaio (LMS) [2^ in comparison 
to runs with the traditional rate treatment. Despite of 
the visible variations with the rate treatment, however, 
the spread of results for different progenitors does not 
widen and again the core properties seem to converge 
during collapse by a self- regulation of electron captures. 
It is unlikely that this result will change when incoherent 
neutrino scattering off nuclei is included in the models. 
The effects of this process during stellar core collapse 
have not been satisfactorily explored yet. 



There is still considerable uncertainty in the supernova 
simulations due to our incomplete knowledge of the nu- 
clear equation of state. The runs for the different pro- 
genitors as well as the studies with varied electron cap- 
ture rates were all performed with the nuclear equation 
of state (EoS) of Lattimer and Swesty (L&S) [2^, which 
is most widely used in core collapse simulations. It is 
based on a compressible liquid drop model and employs 
a Skyrme force for the nucleon interaction. Our choice of 
the compressibility modulus of bulk nuclear matter was 
180 MeV, and the symmetry energy parameter 29.3 MeV, 
but the differences in the supernova evolution caused by 
other values of the compressibility of this EoS were shown 
to be minor pH^Tj. 

We have recently tested the effects of the nuclear EoS 
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FIG. 2: Same as Fig.^but for the rms neutrino energies as defined in Eq. Q. Note that during core collapse and in particular 
before shock formation only i^e were taken into account in the simulations because the production of and heavy-lepton 
neutrinos is suppressed due to the low entropy and correspondingly high electron and Ue degeneracy. 



by using two available alternative descriptions , 
on the one hand a new relativistic mean field EoS 
("Shen") .')(» with a compressibility of nuclear matter 
of 281 MeV and a symmetry energy of 36.9 MeV, and on 
the other hand an EoS which was constructed by Hartree- 
Fock calculations with a Skyrme force for the nucleon- 
nucleon interaction ( "Wolff" ) 31] and has a compressibil- 
ity of 263 MeV and a symmetry energy of 32.9 MeV. In 
particular the latter EoS allows for a faster deleptoniza- 
tion of the less compact and less opaque postshock layer, 
thus producing a clearly higher burst in comparison 
to the standard case with L&S EoS. Note that all three 
runs were performed with the 15 Mq progenitor model 
sl5a28 of Woosley et al. 1^ and included the use of the 
new LMS rates for electron captures. 

Uncertainties of core collapse simulations due to the 
use of different numerical schemes for hydrodynamics and 
neutrino transport were recently investigated in a com- 
parison of the Oak Ridge-Basel and Garching supernova 



codes [131]. Despite of differences in details, very sat- 
isfactory agreement was found for the overall evolution 
and for the neutrino emission properties. The approxi- 
mative description of general relativity in the Garching 
code produces only insignificant deviations from the fully 
relativistic treatment of the Oak Ridge-Basel code during 
the infall phase and early postbounce evolution including 
the prompt neutronization burst. The burst heights and 
widths agree nicely between both codes in Newtonian as 
well as relativistic runs. 

Among the remaining systematic uncertainties in core 
collapse models with possible consequences for the neu- 
tronization burst is the unsettled question of rotation in 
the progenitor core. While the large asymmetries seen in 
supernova explosions are sometimes claimed to be caused 
by rapid core rotation (e.g., Ref. [s^), recent stellar evo- 
lution models seem to disfavor this possibility because 
they predict that massive stars lose angular momentum 
very efficiently during their evolution. The stellar cores 
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therefore rotate so slowly — rotation periods at the edge 
of the iron core before the onset of gravitational insta- 
bility were determined to be around 100 s — that core 
collapse and bounce remain essentially unaffected by ro- 
tational effects [s^l- Instead of rotation, low- mode con- 
vective instabilities have been discussed in this case as 
potential origin of the observed global anisotropics of su- 
pernovae 34, 35] . For a more reliable determination of 
the conditions in the pre-coUapse cores, however, truely 
multi-dimensional stellar evolution models are desired in- 
stead of the currently used spherical ones that are sup- 
plemented by evolution equations for the lateral averages 
of the angular momentum and magnetic field. 

If rapid rotation of the pre-coUapse iron is still con- 
sidered as a viable possibility, despite of probably valid 
objections based on current stellar evolution models, one 
may wonder about the effects of such rotation on the 
prompt Ve burst. For having a noticeable influence, the 
central spin period before collapse must be significantly 
less than 3-5 s, which decreases during collapse by a fac- 
tor of about 100. Unfortunately, all numerical studies 
of rotational core collapse published so far were done 
with very simplistic or no treatment of neutrino trans- 
port (see, e.g., Ref. [l^ and references therein), and only 
the paper by Fryer and Heger [87| | provides information 
in some detail about the neutrino emission, although the 
grey diffusion scheme used in that work falls much be- 
hind the quality and refinement of the transport dis- 
cussed here for simulations of nonrotating (or sufRciently 
slowly rotating) collapsing stars. Besides a possible de- 
pendence of the neutrino signal from the viewing angle 
(as a consequence of the rotational deformation of the 
core and differences of the sho ck p ropagation and break- 
out between pole and equator |33|), the magnitude of the 
neutronization burst and the mean energies of neutrinos 
emitted during the burst seem to be reduced by rapid 
rotation For conclusive results, however, one has 

to await simulations with a better treatment of neutrino 
transport. 

III. ANALYSIS OF THE NEUTRINO SIGNAL 

A. Neutrino fluxes 

The neutrino flux spectra F^. arriving at the Earth are 
determined by the primary neutrino fluxes F^. as well as 
the neutrino mixing scenario, 

Fu^ = +(l-p)^^°^, (2) 

F,^ = pFl + {l-p)Fl , (3) 
4F,^ = (1 - p)Fl + (1 - p)Fl + {2+p + p)Fl , (4) 

where p (p) is the survival probability of an electron (anti- 
)neutrino after propagation through the SN mantle and 
the interstellar medium. We restrict our analysis to the 
standard case of three active neutrino flavors and negli- 
gible magnetic moments or decays |52j . We assume also 



Scenario Hierarchy sin^ ^is p p 
A Normal > 10"^ cos^ j?© 

B Inverted > 10^^ sin^ i9q 

_C Any < IQ-'^ sin^ '&q cos^ 

TABLE I: Survival probabilities for neutrinos, p, and anti- 
neutrinos, p, for various mixing scenarios in case of the den- 
sity profile of the SN progenitor. Terms of the order i?i3 and 
smaller have been neglected. The solar mixing angle is as- 
sumed to be sin'^ i?© = 0.30; for a recent discussion of allowed 
neutrino oscillation parameters see Ref. [4^. 



that the neutrinos do not cross the Earth before reach- 
ing the detector. The main consequence of Earth matter 
effects on ly^ neutrinos — the key channel to observe the 
neutronization burst — is a regeneration effect in scenario 
B and C, thereby slightly improving the chances to detect 
the Ve burst, while Earth matter effects have no impact 
on the signal in scenario A. Therefore, Earth matter ef- 
fects increase the differences between scenario A and B/C 
and it is conservative to neglect them in our analysis. 

The probabilities p and p are basically determined by 
the flavor conversions that take place in the resonance 
layers, where pj-cs « TONAmf cos2i9/(2V2G'fFo£'). Here 
Amf and "i? are the relevant mass difference and mixing 
angle of the neutrinos, ttin is the nucleon mass, Gp the 
Fermi constant and Yc the electron fraction. In contrast 
to the solar case, SN neutrinos must pass through two res- 
onance layers: the H-resonance layer at ^ 10^ g/cnr^ 
corresponds to Am^^j^-^ , whereas the L- resonance layer at 
PL ~ 10 g/cm'^ corresponds to AtUq This hierar- 

chy of the resonance densities, along with their relatively 
small widths, allows the transitions in the two resonance 
layers to be considered independently 

The neutrino survival probabilities can be character- 
ized by the degree of adiabaticity of the resonances tra- 
versed, which are directly connected to the neutrino mix- 
ing scheme. In particular, the L-resonance is always adi- 
abatic and appears only in the neutrino channel, whereas 
the adiabaticity of the H-resonance depends on the value 
of t9i3, and the resonance appears in the neutrino or anti- 
neutrino channel for a normal or inverted mass hierarchy, 
respectively. Tabled shows the survival probabilities for 
electron neutrinos, p, and anti-neutrinos, p, in various 
mixing scenarios for the static density profile of the pro- 
genitor. Using this profile is appropriate, because we are 
only interested in the neutrino propagation during the 
first milliseconds after core-bounce when the shock wave 
has not reached the H-resonance yet (9| . For intermediate 
values of -di^, i.e. 10^^ < sin^iJis < 10^"^, the survival 
probabilities are no longer constant but depend on the 
neutrino energy as well as on the details of the density 
profile of the SN. 

For large values of ??i3, sin^ i?i3 > 10^'^, the H- 
resonance is adiabatic. In the case of a normal mass 
hierarchy, scenario A, the resonance takes place in the 
neutrino channel. Then practically all Ve which are ini- 
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tially created as 1/3 leave the SN also as v^. When they 
reach the detector, they have only a small v^, admixture, 
(i/ajfe) — sin??i3. Taking into account the experimental 
constraints on i?i3, sin^ iJia < 0.047 at 3cr C.L. one 
obtains = sin^ i?i3 < 0.047. This corresponds to an al- 
most complete interchange of the and spectra. For 
an inverted mass hierarchy, scenario B, the resonance oc- 
curs in the anti-neutrino channel, thus interchanging now 
almost completely the 9^ and Dx spectra. In contrast, the 
H-resonance is strongly non-adiabatic for small values of 
i?i3, sin^iJig < 10-^ and for any mass hierarchy (sce- 
nario C), and hence it is ineffective. In both the scenarios 
B and C, the primary v^, leave the star as 1^2 with a large 
Uf, admixture at the detector, (1^2 1 ^^e) = sin-i^Q, leading 
to p^ '^ =sin2|?0 =0.30 [43. 

Let us assume for simplicity that during the neutron- 
ization bursts only i/^ neutrinos are emitted. Then the 
neutrino fluxes arriving at the Earth are in scenario A 

< = + «0, (5) 

2F^^ = (1 - p^)F°^ + (I + p^)FS^ « FO , (6) 

(7) 

and in scenario B or C, 

fZ'"" = + (1 - « sin^ , (8) 

« cos^^qF^^. (9) 

Hence a detector able to identify i/g events will observe 
a peak in the Ve channel in the cases B and C, while 
the peak would be completely absent in case A. The sig- 
nature is similar, although less dramatic, for a detector 
observing elastic scattering events. In this case not only 
but also contribute to the signal through neutral- 
current reactions. But since the cross section for elastic 
scattering on electrons is larger for than for neu- 
trinos, the event number during the neutronization burst 
even for such a non-ideal detector is much larger in the 
scenarios B and C than in A. 

Finally, we note from Fig. ^ that the emission of other 
flavors than i/g becomes important already during the 
end of the neutronization burst, washing out the big dif- 
ferences expected in the naive picture above. In the next 
subsection, we discuss in detail how the neutronization 
burst can be identified. 



B. Detection of the neutronization burst 

Theoretically, the identification of the neutronization 
burst is cleanest with a detector using the charged- 
current absorption of neutrinos. Examples of such 
detectors are heavy water detectors like SNO us- 
ing + d ^ e~ + p + p, or liquid argon detectors like 
ICARUS in using i^^ + '^"Ar + ^^K*. 

The simplest possible observable to identify the neu- 
tronization burst is the total number TV of events 



within an arbitrary fixed period tmax after the onset of 
the neutrino signal. For instance, Ref. I3| calculated 
the expected number of I'e + ''''Ar — > e~ -I- ^^K* events 
in a 70-kton liquid argon detector for tmax — 240 ms, 
where t = 200 ms corresponds to the core bounce, as- 
suming one specific astrophysical model and as SN dis- 
tance d = 10 kpc. They found iV = 41 events in scenario 
A, compared to = 86 in scenarios B and C. From the 
discussion in Sec. In] it is clear that the uncertainty in N 
coming from SN models is rather small. We will quan- 
tify this uncertainty later in Sec. IIII C 21 and use here 
10% as estimate for the systematic uncertainties due to 
the SN models. Combining these systematic uncertain- 
ties and the statistical fluctuations in quadrature leads 
to A^ = 41 ± 8 for scenario A and iV = 86 ± 13 for B 
and C. Thus one could conclude that a differentiation 
between the scenarios A and B/C on the 2 a level is pos- 
sible with a liquid argon detector using the total num- 
ber of events. However, we have still neglected another 
important source of uncertainty for N: The distance to 
stars in our Galaxy is typically known only with 25% 
accuracy The measurement of the SN lightcurve 

will allow a determination of its distance with an error 
of « 5-10% Hi]. However, the probability that the SN 
is obscured by dust is as high as ^75%. Without an es- 
timate for the SN distance, the total number of events 
observed cannot be connected to the SN luminosity and 
is thus not a useful observable. Instead, we exploit in 
the following the time structure of the detected neutrino 
signal as signature for the neutronization burst. 

Since the Ve burst lasts only about 25 ms, cf. Fig.^ the 
event number in current and proposed charged-current 
detectors is not high enough to allow for a detailed time 
analysis. Therefore, we will concentrate in the follow- 
ing on the case of a megaton water Cherenkov detector, 
proposed e.g. to be build in Japan or in the United 
States . A drawback of this choice is that this detector 
type does not have a clean signature for the Ve channel. 
Instead, one has to consider the Ve elastic scattering on 
electrons, i^e + e~ —^ Vg + e~ . This reaction is basically 
affected by three different kinds of backgrounds: inverse 
beta decay reactions De + P n + , reactions on oxy- 
gen, and the elastic scattering of other neutrino flavors 
on electrons. In Fig. |31 we show the distribution of the 
reconstructed energies and directions with respect to the 
vector SN-Earth made of all events in a water Cherenkov 
detector for t < t^ax — 18 ms, where thnunr.p. = ms. 
The events are simulated following Ref. |49| . 



The dominant source of background events are in- 
verse beta decay reactions. These events are almost 
isotropically distributed, while their energy distribution 
reflects the neutrino energy spectrum. In contrast, elastic 
scattering events are concentrated in the forward direc- 
tion and at rather low energies. Therefore a cut with 
E(MeV) X angle(deg) < 500 as shown in Fig. |21 by a 
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FIG. 3: Energy and angular distribution of all events with 
t < tmax ~ 18 ms, assuming tbounco = ms at bounce, in 
a megaton water Cherenkov detector for a SN with a pro- 
genitor mass of 15 Mq at 5 kpc, and scenario C. The dif- 
ferent reaction channels shown are elastic scattering on elec- 
trons (red), inverse beta decay (green), CC events on oxygen 
(black) and NC events on oxygen (blue); also shown is the 
cut _E(MeV) X angle (degrees) = 500. 
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FIG. 4: Number of events per time bin from the elastic scat- 
tering on electrons (black) , inverse beta decay assuming 90% 
tagging efficiency of Gadolinium (red), and reactions on oxy- 
gen (blue) in a megaton water Cherenkov detector for a SN 
with a progenitor mass of 15 Mq at 10 kpc, and case C. Solid 
lines stand for the number of events without cut and dashed 
lines with the cut (_E/MeV) x (angle/ degrees) = 500. 



solid black line substantially reduces the number of back- 
ground events while keeping most of elastic scattering 
events. This background is further reduced by using 
in addition Gadolinium to tag_the neutrons produced 
in the inverse beta reaction 50]. Reactions on oxygen 
have a large reaction threshold, i?th > 15 MeV, and are 
therefore not numerous. Moreover, the charged-current 
events on oxygen (black dots) have an angular distribu- 
tion peaked slightly backwards, and thus the chosen cut 
eliminates these events efficiently, too. 

As illustration for the efficiency of the background sup- 
pression we show in Fig. 0] the expected signal without 
(solid lines) and with cuts (dashed lines). In both cases 
we assumed an efficiency of 90% for the Gd tagging of the 
inverse beta decay reactions. While the number of elas- 
tic scattering events is practically unchanged, the back- 
ground of inverse beta decays and reactions on oxygen 
is significantly reduced. Therefore, we will consider only 
the elastic scattering reactions on electrons in our discus- 
sion below. The sample of elastic scattering events still 
contains the irreducible background of scattering on elec- 
trons of other neutrinos than I'e, but we will show that it 
is possible to disentangle the scenarios with and without 
peak even in the presence of this background. 



C. Results 



In this subsection, we examine in detail how the 
prompt neutronization burst from a future galactic SN 
appears in a water Cherenkov detector. For this purpose, 
we have generated elastic scattering events of neutrinos 
on electrons using Monte Carlo simulations as described 
in Ref. 49]. If not otherwise stated, we have assumed 
a megaton detector with energy threshold i?th = 5 MeV 
and 10 kpc as the distance to the SN. If the energy thresh- 
old could be lowered to 3 MeV, then the event number 
would typically increase by 20%. The neutrino luminosi- 
ties and energy spectra are based on the SN models de- 
scribed in Sec. II. In order to follow the time evolution 
of the signal, we have considered the time interval from 
t — 5 ms until t = 18 ms postbounce and divided the 
interval into five bins. As far as the neutrino mixing sce- 
nario is concerned, we will compare only the cases A and 
C. The first reason for this choice is that the differences 
between B and C arising in the anti-neutrino channel, 
7^ I are always smaller than their differences to case 
A. Secondly, the mixing scenario B can be confirmed or 
ruled out by the modulations in the spectrum induced 
by shock waves in the SN envelope or by Earth matter 
effects, respectively. However, these modulations are not 
helpful to distinguish the cases A and C, since = p'^ . 
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FIG. 5; Number of events per time bin from the elastic scat- 
tering on electrons: total number (black), i/e e (red) and e 
(blue) , for the mixing scenarios C (solid lines) and A (dashed 
lines) . A megaton water Cherenkov and a SN with progenitor 
mass 15 Mq at 10 kpc is assumed. 



1. Dependence on the neutrino mixing scenario 

To understand better the results of our Monte Carlo 
simulations, we first discuss qualitatively the behavior of 
the expected signal in the presence of the irreducible i^^ 
background. The total number of events can be decom- 
posed as 

N{t) = N,^ it) + N,^ it) + N,^ (t) + N,^ (t) , (10) 

where N^. (t) represents the number of events produced in 
the reaction i'i + e~ —^ Vi + e~ . Since we are interested in 
the differences between scenario A and C, and = P^, 
we will concentrate on the number of events produced 
by i^eC" and Vxe~ scatterings. Taking into account the 
linear dependence of the cross sections on the neutrino 
energy, cr^e(-E) « c^e-^j ^'^^ Eqs. we can write 

N,^ (t) cx [pL,^ (t) + (1 - p)L,^ {t)] , (11) 
NuAt) « <J(l-p)L,^(t) + (l+p)L,Jt)] ,(12) 

where crj^ ^ ~ 6cr1 ^. We parameterize our simulation 
results for the different neutrino luminosities by L^^ (t) — 
e{t)Lv^{t), where e{t) is zero until core bounce {tph = 
ms), increases afterwards, and reaches e{t) « 0.5 at tpb — 
18 ms. 



We now examine the difference in the total number of 
events between scenario A and C. We define the following 
ratio. 



_N^At)+N^At) 1 + 7 sit) 



In the first three bins e{t) « and, therefore, all events 
are generated from F° . In scenario C, 30% of the orig- 
inal Ve remain as Ve whereas 70% are converted into 
Vx- In case A, almost all Ve arrive at the Earth as 
h'x- Since the cross section is smaller for i'xS~ than for 
h'eS" , fewer events are expected in scenario A. This dif- 
ference, though, is not constant but changes with s{t), 
from R^'^{t) w 0.4 right before the bounce until 0.9 at 
tph — 18 ms. This evolution can be clearly observed in 
the behavior of the solid (case C) and dashed (case A) 
black lines in Fig. [S] showing the most important con- 
tributions to the neutronization burst signal, t'eC^ (red) 
and VxG' (blue), and the total iye~ (black). 

Next, we study the different contributions from i^e and 
i/x scatterings to the signal. In particular, from Eqs. Ijlll - 
I12|l we define the ratio 



Rex{t) 



6 



p+(l-p) e{t) 
\-p+{l+p) e{t) 



(14) 



iVC (t) + iVC (t) 2.5 + 5.5 e{t) 



(13) 



In case C, R^xi^) > 1 during the whole burst. Therefore, 
i'eS~ scatterings generate always more events than i'xS~ 
reactions, as can be seen in Fig.[Sl This is very different in 
case A: The value of the ratio, R^xit) ~ 6e{t)/{l + e{t)), 
strongly depends on time. In the first two bins, e{t) « 0, 
and thus Rfxi^) ~ 0- Therefore, practically all events 
are generated by Vx^ - At ipb > 11 ms, however, e(t) 
has increased so much that i?^ > 1. In Fig. [3 we can 
see how in the first three bins the dashed blue lines are 
above the red ones, but in the fourth bin both lines cross 
and events from Vee~ scatterings become more impor- 
tant. This feature will play a key role in disentangling 
both mixing scenarios. 

Finally, we discuss the time evolution of the signal. 
This evolution results from an interplay of the time de- 
pendence of (t) and L^^ (t) . Whereas (t) shows a 
clear peak structure around tpb ~ 7 ms, Li,^ (t) steadily 
increases after the bounce. In order to discuss which 
time behavior dominates, we consider separately the two 
channels N^^{t) and N^^{t). From Eqs. (|11II12|) . we can 
estimate the ratio between events connected to L^^{t) 
and to L^^(t). In scenario C, this ratio is smaller than 
one for both channels N,y^ (t) and iV^^ (t) . Thus the com- 
ponent generated by L^^t), and therefore producing a 
clear peak, dominates over L^^{t). This is reflected in 
Fig. O in the peak structure of both the sohd red and 
blue lines. In scenario A, the ratio is smaller than one 
in N,y^ (t) , but it is larger that one in N^^ (t) . There- 
fore, Li^^t) dominates over L^^{t), if i?^(t) < 1. As a 
consequence, the net result in scenario C is an enhanced 
peak, whereas in scenario A the time structure is more 
complicated. In the first three bins, N,y^ (t) dominates 
over N^^t) producing a peak like in case C, but with 
fewer events. In the fourth and fifth bin, however, N^^ (t) 
becomes larger and there is a cancellation between the 
negative slope of N^^t) and the positive of N^^t). The 
final result is an almost flat structure in case A (dashed 
black lines in Fig. O , in contrast to the clear decrease of 
events in the last two bins in case C. 
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FIG. 6: Distribution of fit values a and b for different SN 
distances, D = 2, 5 and fO kpc, and a 15 Mq SN progenitor. 



In order to illustrate the difference in the slope after 
the third bin we have simulated the neutrino signal from 
10000 SN explosions for a progenitor with 15 Mq. Then 
we have fitted the last three bins in each case by a straight 
line, y = at /ms + b. In Fig. |^ we show the normalized 
distribution of the fit values of a and b for scenarios A and 
C and three different distances, D — 2,5 and 10 kpc. The 
slope, a, is centered at almost ms^^ in mixing scenario 
A, whereas in C the center lies at about —3 ms^^. 



2. Dependence on the progenitor mass and SN models 



FIG. 7: Number of events per time bin for different reactions 
in a megaton water Cherenkov detector for a SN at 10 kpc 
for cases A (dashed lines) and C (solid lines) and for different 
progenitor masses: 13 Mq (nl3) in red, 15 Mq (sl5s7b2) in 
magenta, and 25 Mq (s25a28) in blue. Statistical errors are 
also shown for the 15 Mq case. 

progenitor mass from 11 AIq to 25 Mq. This variation 
is at least a factor two smaller than the statistical fluc- 
tuations. However, after the third bin the predictions for 
different progenitor models vary more strongly, reaching 
10% or the size of the statistical fluctuations. At this 
time, neutrinos other than v^. start to play an important 
role and the predictions become more model-dependent, 
cf. Fig. n These differences, though, do not affect the 
observation of the neutronization peak, and therefore do 
not spoil the possibility to disentangle the scenarios A 
and C. 



We have discussed how the presence or absence of a 
peak in the neutrino signal during the neutronization 
burst is related to the neutrino mixing scenario. In this 
subsection, we analyze the robustness of this connection, 
studying its dependence on several SN parameters. 

If the SN is optically obscured, then the progenitor 
cannot be identified. Thus we have to ensure that the 
neutrino signal during the Ve burst depends only weakly 
on the progenitor model. We have analyzed the ex- 
pected neutrino signal for various progenitors with main 
sequence masses 11.2 Mq (si 1.2), 13 Mq (nl3), 15 Mq 
(sl5s7b2)M5 Mq (sl5a28), 20 Mq (s20) and 25 Mq 
(s25a28) Ji3|. In Fig. we show for illustration the re- 
sults for two extreme progenitors, nl3 and s25a28, as well 
as for an intermediate case, sl5s7b2. 

We find that for all progenitor masses the peak struc- 
ture described in the last subsection can be clearly ob- 
served in case C but not in case A. In the first three 
bins the differences between different progenitor masses 
are smaller than the statistical ffuctuations. The average 
event number per bin varies less than 6% changing the 



Another source of uncertainty in our predictions of 
the neutrino signal may be the incomplete or inaccurate 
treatment of some of the weak interaction rates, in par- 
ticular with nuclei, used in the SN simulations. As an 
example for the changes in the neutrino signal induced 
by an improved treatment of these rates we consider the 
effect of electron captures by nuclei in the core. We com- 
pare the signal expected for the progenitor models sl5a28 
and s25a28, see Fig. |S1 with and without including the 
electron capture rates on nuclei with mass numbers larger 
than 65 as suggested in Ref. ^ (LMS). In both sl5a28 
and s25a28 SN models we observe the main differences in 
the central bins. This change is directly connected to the 
behavior of the luminosities: Electron capture on heavy 
nuclei broadens the peak in L,^^ (t) and reduces somewhat 
{t). However, the changes are always smaller than the 
statistical fluctuations and therefore do not affect the ob- 
servation of the neutronization peak. 
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FIG. 8: Same as Fig.lZlfor mixing scenarios A (dashed lines) 
and C (solid lines) for a SN at 10 kpc and a progenitor mass 
of 15 Mq (sl5a28). In blue we show the SN model with the 
electron capture rates on heavy nuclei adopted from Ref. |23l 
(LMS), and in red with the traditional description. Statistical 
errors are also shown for the (LMS) models. 
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FIG. 9: Same as Fig.Elfor mixing scenarios A (dashed lines) 
and C (solid lines) for a SN at 10 kpc and a progenitor mass of 
15 Mq (sl5a28). We show the SN models computed with the 
equations of state "L&S" (red), "Shen" (blue), and "Wolff" 
(magenta), respectively. Statistical errors are displayed for 
the "Shen" model. 



For the SN progenitor sl5a28 we have also investigated 
the influence of the nuclear equation of state (EoS) on 
the evolution of the luminosity during the neutronization 
burst. We have considered three different models for the 
EoS, dubbed "L&S" [H, "Shen" |30|], and "Wolff" 
In Fig. 1^1 we show the predicted numbers of events de- 
pending on the employed EoS. The main change is the 
total luminosity in the peak. In case C, the change in 
the luminosities results only in a rescaling of the total 
number of events in each bin. In case A, though, the 
increase of L^,^ leads to an enhancement of the + e~ 
contribution, whereas the Ve + e~ events are not strongly 
affected. As can be seen in Fig.|31 this may slightly mod- 
ify the time-structure and may lead to a small peak in 
the third bin. However, these changes are again much 
smaller than the size of the statistical fluctuations. 



In summary, we have found that the changes in the 
predicted event numbers for different progenitor models, 
improvements in the interaction rates or changes in the 
EoS are small compared to the size of the statistical fluc- 
tuations. However, up to now we have always varied 
only one parameter, e.g. the progenitor mass, and kept 
all others fixed. In order to quantify the probability to 
disentangle the neutrino mixing scenarios A and C, we 
should in principle first quantify the uncertainties in all 
input parameters of the SN model and then derive the 
resulting uncertainties in the neutrino fluxes. Already 
the first step is impractical. Therefore, we use the fol- 
lowing method: We take all the SN models considered in 
the previous sections and calculate for each the expected 



number of events for cases A, N^^, and C, N^,^, where 
a refers to a SN model. Then we compute the proba- 
bility that scenario A and C can be distinguished for all 
possible pairs a,/? using a analysis, 

5 (/VA._/VC.)2 
Xac(Q^./^) /VA 4, /VC ■ 

i=l ^^Q,j + ^^/3,i 

In Fig. ^] we plotted with solid lines the distribution of 
this probability using all possible combinations of the SN 
models introduced previously. In the worst combination, 
scenario A and C can be distinguished only at the 2a 
level {p = 4.5%). However, in most cases the probabil- 
ity to misidentify scenarios A and C is smaller than 3ct 
[p = 0.27%). We show in Fig. [TOI also the case that the 
systematic uncertainty of the SN models is dominated 
by the unknown progenitor star (dashed line). Then it 
is possible to disentangle the scenarios A and C with a 
confidence level better than Scr for most combinations. 



3. Intermediate values of-dis and scenario B 

For completeness, we analyze in this subsection the 
dependence of the neutrino signal on intermediate val- 
ues of i?i3 in scenario A, and the case of scenario B. For 
this purpose we fix a SN progenitor model, sl5a28, and 
consider first mixing scenarios that lie between A and 
C, i.e. a normal mass hierarchy and sin^ t?i3 = 10~^, 
5 X IQ-^ 2 X 10"^ and 10-^ see Fig. [TT] The first and 
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FIG. 10; Distribution of the probability to misidentify cases 
A and C for all SN models considered in the previous sec- 
tions combined (solid lines), and only for those with different 
progenitor mass (dashed lines). 



the last value correspond to scenario A and C, respec- 
tively. Whereas p{E) is independent on i9i3, the average 
value of p{E) grows continuously from in mixing sce- 
nario A to sin^ 'i?0 ~ 0.3 in scenario C. This behavior of 
the survival probability results in a continuous increase 
of the peak. For instance, the probability to disentangle 
the case of intermediate iJis from scenario A becomes 
smaller than 2(t for sin^ -di^ > 10^^ and the SN progen- 
itor sl5a28. Therefore, the detection of a peak not only 
excludes case A (sin^i?i3 > 10^'^), but also most of the 
intermediate range of z?i3, sin^i^ia > 10"^ 

Finally, we show also an example of case B, inverted 
mass hierarchy and large t^ia, sin^ -i^ia — 10^^. Since the 
survival probability for anti-neutrinos is different for B 
and C, the total number of events also differs slightly. In 
particular, the number of events in case B is expected 
to be larger, as it can be seen in Fig. ^2 The reason is 
twofold: First, Lp^ > Lp^, and second a^^e > fp^e- Since 
Ki 0, more Dx are converted into De than in case C, 
and the larger e~ cross section implies more events. 
Hence, we conclude that scenario B is slightly easier to 
distinguish from C than A from C. 



IV. DISTANCE DETERMINATION 

After the neutrino mixing scenario has been deter- 
mined, the small spread in the predicted total number 
of events N suggests that is a useful estimator for the 
SN distance D. Since D cx N^^"^, the relative uncertainty 



FIG. 11; Same as Fig. 13 for a SN at 10 kpc and a progen- 
itor mass of 15 Mq (sl5a28), for different neutrino param- 
eters. We show with a dashed line case B, inverted mass 
hierarchy, and with solid lines a normal mass hierarchy for 
different values of sin^i?i3. Statistical errors are shown for 
sin^tJis = 2 X 10-^ 



SN/Nq translates into 



SD=^6N. 
2No 



(16) 



For each SN model a, the probability distribution Pa{N) 
of the observed number of events N is Gaussian with 
CTstat = \/ (Na). To estimate the "systematic" uncer- 
tainty (Tsyst due to differences in the SN models, we use 
the following recipe: We calculate for different models 
the individual expectation value (Na) and use then their 
variance as systematic uncertainty, i.e. 



1 



sys 



1 ^ 



(17) 



where m stands for the total number of SN models con- 
sidered, and (N) = (l/m)X;„(^a)- 

The two most important sources of systematic errors 
are the nuclear equation of state, af°^ and, if the SN pro- 
genitor is not identified, its mass, cr^s^^- In order to quan- 
tify the influence of different progenitor types we have 
considered six different SN models (sll.2, nl3, sl5s7b2, 
sl5a28, s20, and s25a28) discussed in Sec. II. For the time 
interval considered here, from t — —5 ms until t = 18 ms 
post bounce, we obtain as systematic error <j"ys^^ ~ 7%. 
In the case of a SN located at 10 kpc we expect as aver- 
age value (N) = 112 and o'sytt'' w 8 as systematic error 
for scenario A, compared to {N) = 176 and cr^ft*' ~ 12 
for scenario C. In Fig.^J we show the probability distri- 
bution Pa (N) of the observed number of events N until 
t = 18 ms post bounce, choosing as examples three SN 
models: sl5s7b2, and the extreme cases nl3 and s25a28. 
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FIG. 12; Distribution pa{N) of the observed total number of 
events A*" for 20000 simulated SNe with progenitors 

: nl3, sl5s7b2, and s25a28, and for the neutrino mixing 
scenarios A and C. A distance to the SN of 10 kpc has been 
assumed. 



FIG. 13: Same as Fig.[T2lfor the SN progenitor sl5a28 and 
three different EoS: L&S, Shen and Wolff. 



The obtained range [(TV) - cr^f' : (N) + a^f""] is also 
shown by arrows for scenarios A and C. 



The second source of systematic uncertainty analyzed 
is the nuclear equation of state. As an illustration of its 
effect on the determination of the SN distance we have 
considered the SN progenitor sl5a28 using the three dif- 
ferent EoS defined in Sec. II: L&S, Shen and Wolff. While 
the energy emitted in v^. and until i = 18 ms post 
bounce in the models L&S and Shen differs by less than 
3%, the energy released in the Wolff model differs by 
more than 16%. This difference depends on the flavor: 
in the Wolff model more I'e's but fewer v^^s are emitted 
than in the other models. In scenario C the signal is 
dominated by L^,^ , and thus the number of events in the 
Wolff model is much larger than in the other two mod- 
els. However, in scenario A the contribution from 
is of the same order of magnitude as that of L^^ , and 
since L,y^ is smaller in the Wolff model, this compen- 
sates the increase in the number of events due to L^^. 
This leads to much less variation in the total number 
of events in case A than in C. For the same time inter- 
val we have obtained a systematic error af"^ « 2% and 
6% in cases A and C, respectively. If a SN is located 
at 10 kpc we expect as average value (N) ~ 106 and 
af°^ ~ 2 as systematic error for scenario A, compared to 
(A^) = 189 and cr|?J « 12 for scenario C. In Fig. [T21 we 
show the probability distribution pa {N) of the observed 
number of events N until t = 18 ms post bounce for 
the three different models considered. The corresponding 
range [{N) — af°^ : (N) + crfy^] is also shown by arrows 
for scenario A and C. 



Assuming that cr^^^^ and 0"^° dominate the total sys- 
tematic uncertainty (Tsys we obtain a total "systematic 
error" of 7% and 9% for scenarios A and C, respectively. 
We can now combine the statistical errors associated with 
the averages (N) and systematical errors again in quadra- 
ture. For the cases discussed in this section this results 
in a relative error of 12% for the total event number, or a 
6% error in the SN distance as predicted from the prompt 
neutronization burst. 

This result represents a significant improvement with 
respect to the typical accuracy of the distance determi- 
nation for stars in our Galaxy, which is of the order of 
25%. In the case that the SN lightcurve can be measured, 
the distance to the star could be determined with a sim- 
ilar accuracy, 5-10%. However, it is likely that the next 
Galactic SN is obscured by interstellar dust: about three 
out of four SNe are estimated to be optically obscured by 
dust In this case the measurement of the neutrino 
signal from the neutronization burst would provide the 
only way to determine the distance to the star. 



V. SUMMARY 

We have performed simulations of SN explosions for a 
variety of progenitor models and differences in the un- 
derlying microphysics used. The resulting changes in 
the neutrino signal from the neutronization burst were 
always small compared to the statistical fluctuations ex- 
pected for a megaton water Cherenkov detector. We have 
argued that the total number of events is not a useful ob- 
servable, if the SN is optically obscured by dust. Instead, 
we propose to use the time structure of the neutrino sig- 
nal as discriminator for the neutrino mixing scenario: the 
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observation of a peak in the first milliseconds of the neu- 
trino signal would rule out the normal mass hierarchy 
with "large" iJis (case A), cf. Figs. 6-9. Thereby, the 
degeneracy between scenarios A and C, which both pro- 
duce the same signal during the accretion and cooling 
phases, can be broken. 

After the neutrino mixing scenario has been estab- 
lished, the small uncertainty in the predicted value of the 
total event rate during the neutronization burst phase 
makes a determination of the SN distance feasible inde- 
pendent from optical observations. Provided that the 
predictions of current progenitor models and state-of- 
the-art simulations do not miss important physics, we 
estimate that the SN distance can be measured with an 
accuracy of ~6%. Since it is likely that a Galactic SN 
is optically obscured by interstellar dust and no other 
distance determination can be used, this new method re- 



lying only on neutrinos seems very promising. 
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